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Résumé 

Plusieurs indices suggèrent qu'un grand nombre d'astéroïdes seraient des empilements de gravats, liés uniquement 
par gravité et quasiment dépourvus de cohésion interne. Leur formation serait due à la réaccumulation de fragments 
issus de la destruction antérieure d'un corps parent. La forme de ces objets, que l'on observe, pourrait ainsi être liée 
à ce processus de réaccumulation gravitationnelle. Toutefois, à l'heure actuelle, on ignore si les formes observées 
des astéroïdes sont le résultat de la seule réaccumulation ou d'événements ultérieurs entraînant un réajustement 
de la forme du corps. 

Dans cet article, nous nous interrogeons sur l'origine des formes des astéroïdes. A l'aide d'une étude numérique 
de la réaccumulation, nous étudions les formes d'équilibres possibles correspondants aux modèles théoriques. Les 
résultats montrent, pour la première fois, que seulement une catégorie précise de formes (sphéroïdes aplatis) 
est apparemment crée via ce mécanisme. Ces résultats pourraient ainsi fournir d'intéressantes contraintes sur 
l'évolution des formes d'astéroïdes, notamment pour ceux possédant un, ou plusieurs, satellites. 

Pour citer cet article : J.-F. Consigli et ai, C. R. Physique ... (2006). 

Abstract 

Asteroid shapes and satellites: rôle of gravitational reaccumulation. FoUowing current évidences, it is 
widely accepted that many asteroids would be "gravitational aggregates", i.e. bodies lacking internai cohésion. 
They could mainly be originated during the catastrophic disruption of some parent bodies, through the grav- 
itational reaccumulation of the resulting fragments. The same events produced the dynamical familles that we 
observe. In this work we address the problem of the origin of shapes of gravitational aggregates, that could contain 
signatures of their origin. We use a N-body code to simulate the coUapse of a cloud of fragments, with a variety of 
initial velocity distributions and total angular momentum. The fragments are treated as inhelastic sphères, that 
rapidly accumulate to form rotating aggregates. The resulting shapes and rotational properties are compared with 
theoretical prédictions. The results show that only a précise category of shapes (flattened spheroids) are created 
via this mechanism. This may provide interesting constraints on the évolution of asteroid shapes, in particular for 
those with one or more satellites. 

To cite this article: J.-F. Consigli et al., C. R. Physique ... (2006). 



Key words: Solar System: Asteroids: Satellites; Collisions 
Mots-clés : Système Solaire; Astéroïdes; Satellites; Collisions 



Prcprint suhmittccl to Elsevier Science 



2 février 2008 



1. Introduction 

Depuis leur formation, tous les astéroïdes évoluent principalement sous l'influence des multiples perturbations 
planétaires et des chocs qu'ils subissent entre eux. Plusieurs études [1],[2] montrent qu'hormis les plus grands 
d'entre eux (Cérès, Pallas, Vesta, Juno...), les astéroïdes seraient le résultat de collisions mutuelles. 
Parmi ces collisions, les plus violentes d'entre-elles, appelées collisions catastrophiques, sont responsables de la 
destruction totale du corps impacté et génèrent ainsi des centaines de milliers de fragments dont il demeure, encore 
aujourd'hui, des traces observationnelles. Le témoignage de cette intense évolution coUisionnelle est apporté par 
des groupes d'objets partageant des compositions et des propriétés orbitales similaires. L'existence de ces groupes 
d'objets, les familles dynamiques d'astéroïdes, est désormais clairement établie [2]. A ce jour, une vingtaine de ces 
familles a ainsi été découverte dans le Système Solaire. 

Contrairement à ce que l'on pensait jusqu'alors, on sait depuis peu [5] que la formation d'une famille d'astéroïdes 
ne peut s'expliquer uniquement par la fragmentation d'un corps "parent" . Dès qu'un objet dépasse quelques cen- 
taines de mètres, sa destruction signifie qu'il se fragmente lors de la traversée de l'onde de choc, mais également 
que les fragments produits s'échappent et interagissent mutuellement. Certains fragments peuvent se réaccumuler 
sous leur propre attraction pour former des agrégats qui constitueraient ainsi les membres des familles d'astéroïdes 
observées aujourd'hui [5]. 

Outre la formation de ces familles, ces collisions peuvent également entraîner la formation de satellites [3], [6] 
en orbites stables [7], [9] autour d'astéroïdes. Pratiquement dépourvus de toutes forces de cohésion interne, ces 
corps réaccumulés seraient des agrégats de fragments liés uniquement par la gravité mutuelle s'exercant entre les 
différents fragments constituant l'agrégat lui-même. 

Parallèlement à ces travaux, des preuves indirectes de l'existence de ces agrégats ont été obtenues. On peut citer 
à ce propos les études de courbes de lumière [10], [11], ou les observations de l'astéroïde Mathilde par la sonde 
NEAR. Mathilde possède en effet une faible densité, p =1.3 g cm~'^, soit environ 3 fois inférieure aux astéroïdes 
de composition similaire [16]. Un tel constat laisse penser que Mathilde serait en réalité un corps réaccumulé, de 
grande porosité (pourcentage de vide), formé de débris issus d'une collision antérieure. 

L'objet de ce travail est précisément d'examiner les caractéristiques des astéroïdes se formant par réaccumulation 
gravitationnelle lors de l'effondrement d'un nuage de particules. La question centrale que l'on se pose est de 
savoir quels liens il existe entre la forme des astéroïdes, la présence éventuelle de satellites et le processus de 
réaccumulation lui-même. Afin d'apporter des éléments de réponses, on réalise une étude numérique ainsi qu'une 
comparaison avec des modèles de formes d'équilibres déjà existants. L'approche utilisée au cours de ce travail 
se caractérise par son aspect novateur. Une étude précédente [29] prenait en effet comme point de départ des 
agrégats déjà formés et s'intéressait à la façon dont la forme de ces corps évoluait et se réajustait en fonction de 
leur vitesse de rotation. La démarche utilisée ici est différente : on souhaite comprendre quelles sont les formes 
"possibles" (accessibles) via le processus de réaccumulation, et ce, sans faire appel à un corps déjà formés, mais 
en considérant uniquement l'effondrement d'une distribution dispersée de particules, censée représenter un nuage 
de débris après une collision catastrophique. 

Nous évoquerons dans la partie 2 les différents modèles de formes d'équilibres. Dans la partie 3, nous discuterons 
de la méthode numérique utilisée pour mener à bien les différentes simulations. Enfin, en 4, nous analyserons les 
résultats issus de ces mêmes simulations afin d'en tirer des conclusions préliminaires quant à l'origine des formes 
des astéroïdes que l'on observe. 



2. Formes d'équilibre 

Tout ellipsoïde est défini par ses 3 demi-axes a > h > c on par ses rapports d'axes (généralement b/a et c/a). 
Scion certaines études [17], [18] les courbes de lumière observées seraient compatibles avec des formes ellipsoïdales 
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à 3 axes et, plus généralement, avec des formes convexes en rotation uniforme. Récemment, ces résultats ont 
été renforcé pour certains objets grâce à des observations menées avec le télescope spatial Hubble [19], [20]. Ceci 
n'exclut pas le fait que les courbes de lumière puissent contenir, au moins dans quelques cas, beaucoup plus de 
renseignements sur des formes davantage complexes [18]. 

Il y a une vingtaine d'années, certaines études [2], [22] tentèrent de modéliser les astéroïdes ré-accumulés, ou 
^^rubble pile" , par des fluides en équilibre hydrostatique. Cependant, l'approche "fluide", qui consiste à calculer 
la surface équipotentielle pour un fluide en rotation, comporte certaines faiblesses. Si l'on considère les rapports 
d'ajces des astéroïdes observés par photométrie, on s'aperçoit (fig. 3) que peu d'entre-eux possèdent des formes 
d'équilibre de fluides [26], [17]. Cependant, parmi les astéroïdes les plus proches aux formes d'équilibres de fluides, 
on s'est aperçu [27] que les astéroïdes binaires étaient particulièrement nombreux. La réciproque, par contre, n'est 
pas vraie [27] : les astéroïdes binaires connus ne possèdent pas tous des formes d'équilibres de fluides. Dans les 
années 80-90, on s'intéressa davantage aux solutions d'équilibres dans le cas de fluides compressibles, c'est-à-dire 
des fluides dont l'indice polytropique n est tel que n 7^ [24], [25]. Plus récemment, il a été montré [28] qu'un 
faible frottement interne suffisait à rendre compte de quasiment toutes les formes d'astéroïdes observées. 
Nous reviendrons sur ces modèles (2.3) après avoir montré les éléments qui conduisent au calcul des formes 
d'équilibre de fluides (2.2). 



2.1. Généralités et rotation limite 



Tout ellipsoïde triaxial est une surface quadratique dont l'enveloppe est donnée, en coordonnées cartésiennes, 
par : 

Le triplet {a,b,c) représente, respectivement, les demi-axes de l'ellipsoïde : grand, moyen et petit. L'ellipsoïde 
devient une sphère si ses trois demi-axes sont égaux. Si deux des trois demi-axes sont égaux, l'ellipsoïde est appelé 
sphéroïde (ou ellipsoïde biaxial). Selon si c < a ou c > a, il s'agira respectivement d'un obloïde ou d'un proloïde. 
On déflnit ensuite les rapports d'axes b/a et c/a. Sauf mention contraire, l'aplatissement de tout ellipsoïde sera 
défini par e = y^l — {c/a)'^. 

Quasiment dépourvus de cohésion interne, les agrégats sont très sensibles aux forces centrifuges inhérentes à 
la rotation. Dès lors, on peut en effet se demander quelle est la vitesse de rotation limite Qmax à partir de laquelle 
une particule située en périphérie d'un objet sphérique (et plus généralement ellipsoïdal) peut s'échapper sous 
l'effet de la seule force centrifuge. En égalant gravitation et forces centrifuges, on obtient dans le cas d'une sphère 
la relation suivante [29] : 



^y2nGprp V 3 

où prp est la densité du corps considéré. Pour un proloïde, cette relation prend la forme suivante [29] 



— exp ' (e — 1) ' 



^j2-KGprp 



2ë + ln 



1-f ë 



(3) 



avec ici l'aplatissement e est défini par e — ^1 — q| 7^ e. Numériquement, il est parfois utile d'assimiler le proloïde 
à un corps ponctuel de même masse. On aboutit alors à la relation simplifiée [29] : 

— = 7273(1 -ë')-'/'' (4) 

sJï-KGprp 



2.2. Equilibre hydrostatique 



Toute masse fiuide, en rotation rigide, soumise uniquement aux forces gravitationnelles et de pression, suit 
l'équation de l'équilibre hydrostatique : 

p^ = -\7p + pg+^p\7\nx x\'^ + 2puxD, (5) 

on = + (u.V), u est la vitesse au point x, p la pression, g l'accélération de la pesanteur et Q. la vitesse 
angulaire. Dus à la rotation, les deux derniers termes de l'équation sont respectivement la force centrifuge et la 
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force de Coriolis. 

Pour un sphéroïde de Maclaurin (a = b) on peut montrer [23] qu'il existe une relation liant l'aplatissement e du 
sphéroïde à sa vitesse angulaire fl : 

^^2(3-2e^)^arcsin(e)-^ (6) 

De même, il est possible d'exprimer le moment angulaire L en fonction de la vitesse de rotation fi [23] : 

\/3 /a\2 



L = 



avec L = L/{Gm^R) le moment angulaire adimensioné, R = {a^c)^^^ le rayon d'une sphère de même masse m 
que le sphéroïde et fi = Q/{TvGp)^^^ la vitesse de rotation adimensionée. 

Pour un ellipsoïde de Jacobi [a ^ b), on obtient grâce à une algèbre plus compliquée [23] : 

/S.'^ = {a^ + u){b'^ + u){é' + u) (9) 
Là encore, on peut établir une relation liant moment angulaire L et vitesse de rotation fi : 

V3aN_6^^ (10) 

Augmenter la rotation contraint le corps fluide à prendre une forme de plus en plus aplatie. Pour des faibles 
vitesses de rotation, deux formes d'équilibre sont possibles : l'une avec un aplatissement e < 0.8 et telle que a — b 
(sphéroïde; séquence de Maclaurin), l'autre avec la même vitesse de rotation mais un aplatissement plus grand 
(0.8 < e < 1) et telle que a ^ b (ellipsoïde à 3 axes; séquence de Jacobi). Celle-ci possède aussi un plus grand 
moment d'inertie, donc un moment angulaire L plus élevé. Au-dessus de ficrit, aucune forme stable n'existe pour 
un fluide (fig. 1). 

2.3. Modèle élasto-plastique 

Bien que l'hypothèse d'un comportement fluide semble valable pour les corps gazeux du Système solaire, rien 
n'est moins siir à propos des petits corps (astéroïdes) pour lesquels le comportement fluide est loin d'être ga- 
ranti. A priori séduisante, l'idée d'assimiler les astéroïdes à des corps linéairement élastiques comporte quelques 
inconvénients. En effet, les solutions "élastiques" reposent sur une hypothèse forte : l'existence d'un état initial 
sans contraintes dans le passé évolutif du corps considéré [28]. Or, les états actuels des corps du Système Solaire 
sont le résultat d'une évolution mouvementée : collisions, fragmentations, accumulations et réajustements. Ces 
phénomènes sont à l'origine de contraintes résiduelles qui ne peuvent être connues avec précision. 
Adopter une approche élasto-plastique permet de s'affranchir de la connaissance du passé de l'objet. En effet, en 
calculant le chargement maximal accessible pour un corps donné, on peut déterminer des solutions d'équilibres li- 
mites et ce, en particulier, pour des corps composés de fragments. D'après des études géologiques, le comportement 
de certains matériaux granulaires (roches, sols ou graviers) nécessite des théories élasto-plastiques et répond au 
critère de Mohr-Coulomb. Déterminé uniquement par une cohésion Y et un angle de frottement interne </>, ce critère 
est généralement utilisé pour calculer la contrainte maximale qu'un corps peut supporter avant cassure. Etant 
symétrique, le tenseur des contraintes E = [aij] est diagonalisable et admet pour valeurs propres ai > (J2 > o"3, 
appelées contraintes principales. Dans le cas où la cohésion est négligeable, le critère de Mohr-Coulomb s'écrit 
[28] : 

tan(<^) > ^lp=l (11) 

L'angle de frottement est alors limité uniquement par les contraintes principales ai et a^. Typiquement, les 
matériaux granulaires terrestres sont tels que (p ~ 30°. Ainsi, un tas de sable aura pour pente maximale (p ~ 30° 
alors qu'un fluide (0 — 0°) ne formera aucun tas. Il est relativement fréquent d'assimiler un matériau à un solide 
élasto-plastique sans cohésion interne. Couramment employé en mécanique des sols, ce modèle semble approprié 
pour décrire la structure interne des astéroïdes réaccumulés au sein desquels les forces de cohésion sont négligeables 
[28]. 
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3. Méthode numérique 



Afin d'aborder la réaccumulation gravitationnelle d'un système de particules et la formation d'un ou plusieurs 
agrégats de fragments, nous avons étudié l'évolution d'un nuage de particules en imposant une dissipation par im- 
pacts de l'énergie cinétique du système. Cela a été rendu possible par l'utilisation d'un code Af-corps hiérarchique 
s'inspirant fortement de l'algorithme de Barnes-Hut [30], pkdgrav (Parallel K-D tree GRAVity code). Les si- 
mulations ont été effectuées sur une station de calcul parallèle - SIVAM II - à l'Observatoire de la Côte d'Azur. 
Bien qu'effectuées sur une station de calcul parallèle à quatre processeurs, chaque simulation a nécessité, en règle 
générale, plusieurs heures de temps CPU. 

Contrairement aux codes N-corps classiques, pkdgrav a été modifié pour permettre de traiter les collisions entre 
les difi'érentes particules du système. Les particules sont ainsi considérées non pas comme des entités ponctuelles, 
mais comme des sphères de rayons finis. Lorsque deux particules viennent à se rencontrer elles ne fusionnent pas 
mais peuvent rebondir les unes sur les autres conformément aux lois gouvernant les collisions inélastiques entre 
sphères. 

Prenons l'exemple d'une particule P dont les composantes du vecteur vitesse sont et v^. Lorsque la particule 
P entre en collision avec une particule Pl, celle-ci se retrouve, après collision, avec deux composantes de vitesse 
= e||î;|| et vii_ — e±v± où ey et e± sont respectivement les coefficients de restitution parallèle et perpendiculaire 
associés aux collisions inélastiques. 

Dans toutes les simulations, nous avons utilisé une valeur e± = 0.8, ce qui permet une dissipation rapide de 
l'énergie cinétique lors des collisions mutuelles entre particules. Le nombre de particules utilisé {N — 500) s'avère 
être un compromis entre temps de calcul et représentation de la forme du corps final obtenu en fin de simulation. 
[29] a en effet montré dans ses simulations qu'une valeur N > 200 n'introduisait aucun biais important dans la 
représentation (la forme) des objets créés. 

Une fois le corps formé, le code continue à calculer les collisions mutuelles entre particules (même si les sphères 
sont presque en contacts), rendant ainsi la simulation très coiiteuse numériquement. 

3.1. Conditions initiales 

Les conditions initiales consistent en un nuage de masse Mtotai, de N particules de densité pp. Sa taille est 
défini par 3 demi-axes de longueurs dx, dy et d^. 

Pour les simulations numériques présentées ici, nous avons pris Mtotai = 4 x 10^^ kg. Dans le cas idéalisé oii 100% 
de la masse initialement disponible est réaccrétée pour former le corps final, ceci correspond à un rayon des corps 
réaccumulés de Rrp ~ 1 km. Notons toutefois que certaines simulations ont été réalisé avec des densités pp > 1 g 
cm~^ (resp. pp = 2, 3 et 4 g cm"''). 

Un grand nombre de simulations a été réalisé avec dx = dy = dz = 10*m ~ lOTÏrp, c'est-à-dire avec un nuage initial 
10 fois plus grand que le corps final obtenu. Enfin, différents profils de vitesses pour les particules du nuage ont 
été utilisés : un profil uniforme et un profil en 1 /r pour lequel la vitesse des particules diminue vers la périphérie 
du nuage. 

Les différentes conditions initiales sont donc caractérisées par les paramètres suivants : taille et forme du nuage, 
distribution des vitesses et moment angulaire total. 

Une fois initialisé, le nuage de particules est tel que : 

1) les particules qui le composent occupent des positions aléatoires. 

2) les vecteurs vitesse des particules ont un module aléatoire v — rj.Vmax, où r; est une variable aléatoire com- 
prise entre et 1, de moyenne 0.5, et distribuée soit uniformément soit selon une loi en 1/r. On a, de plus, 
< Vmax < k-viib OU k < 1 cst Un paramètre ajustable et vm, est la vitesse de libération du nuage pour une 
particule située à sa périphérie. 

3) Les directions des vitesses sont ajustées aléatoirement par un processus de Monte-Carlo afin d'obtenir la valeur 
de L souhaitée. 

On peut donc considérer les vitesses initiales au sein du nuage comme la somme de deux composantes : une 
composante aléatoire et une autre composante, systématique, qui détermine la rotation d'ensemble du nuage. Ceci 
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permet de sélectionner un "degré de chaoticité" des vitesses initiales, et de reproduire ainsi différentes situations 
rencontrées au cours de simulations de réaccumulation, tels quelques cas typiques de familles d'astéroïdes comme 
"Eunomia" (effondrement rapide) ou "Koronis" (grande chaoticité, effondrement lent) [5]. 



4. Résultats et analyse 

L'agrégat principal est identifié grâce à la proximité géométrique des particules. La vitesse de rotation de 
l'agrégat est obtenue en calculant le moment angulaire total L de l'agrégat, qui lui-même est obtenu en addition- 
nant le moment angulaire de chaque particule prise séparément : L — Lp. Les dimensions de l'agrégat sont 
obtenues en déterminant les demi-axes a, b et c. Pour déterminer a, on considère la distance entre le centre de 
masse de l'agrégat et le bord extérieur de la particule la plus lointaine. Un calcul analogue dans les directions 
perpendiculaires à a fournit successivement une valeur de b et c. Le volume Vrp de l'agrégat {Vrp = | Trafic), permet 
de déterminer une valeur de la densité moyenne prp moyenne. 

L'ensemble des simulations effectuées est présenté dans les différents tableaux (tab. 2 & 3). Les principaux résultats 
sont représentés sur les figures 3 et 4. Ces différentes simulations ont permis d'examiner, en fonction du moment 
angulaire L, les différentes issues possibles du processus de réaccumulation. Nous avons ainsi tenté de vérifier si, 
pour des valeurs élevées du moment angulaire L, il était possible d'observer la formation direct(|3] d'astéroïdes 
binaires via réaccumulation. Certaines simulations décrites ci-après se sont achevées par la formation de systèmes 
binaires, confirmant ainsi la possibilité de former de tels couples par processus de réaccumulation. Les périodes 
de rotation obtenues sont également en accord avec les périodes maximales de rotation définies par les relations 
(2) et (4). 

Au cours des simulations, il est apparu que les paramètres caractéristiques des corps réaccumulés (masse, 
rapports d'axes et vitesse angulaire) convergent très rapidement, et ne permettent pas d'observer une réelle 

"relaxation" graduelle. Le temps dynamique d'effondrement du nuage initial est donné par t^ff = \\J où 
P" ~ 4 . est la densité du nuage initial. L'application numérique donne te// ~ 0.56 jours, ce qui est 

TTa^; Ciy 0,2, 

compatible (en ordre de grandeur) aux 1.45 jours obtenus expérimentalement. 
4.1. Des sphéroïdes aux ellipsoïdes 

Le phénomène de réaccumulation est par essence très chaotique. Il est donc difficile de prévoir a priori la position 
exacte qu'occupera, en fin de simulation, le corps final dans le plan {L, Ù^). 

Pour des faibles valeurs du moment angulaire (Z < 0.2), les simulations montrent que les corps formés sont 
typiquement des sphéroïdes, c'est-à-dire des corps tels que a ~ 6 (fig. 4). De formes sphéroïdales, les corps formés 
se situent au voisinage de la séquence de Maclaurin et possèdent des périodes de rotation P > 10 h. Ces résultats 
tendent à montrer que la phase de réaccumulation post-collision conduit préférentiellement à la formation de corps 
sphéroïdaux à rotation lente ~ 10 h (relativement lente pour des objets d'environ 1 km de rayon). Hormis les cas 
oii l'on assiste à la formation de satellites, ces simulations aboutissent systématiquement à des corps réaccumulés 
à plus de 90 %. Ces corps possèdent ainsi des rayons iî ~ 1 km et des densités de l'ordre de prp — 0.42 g cm~^. 
Des résultats similaires ont été obtenus en augmentant la densité des particules ou, autrement dit, la masse 
du nuage initial. De nouveau, les corps formés en fin de simulation sont proches de sphéroïdes. Comme évoqué 
précédemment, ceci s'explique par le fait que la densité pp des particules ainsi que la masse Mtotai du nuage 
n'infiuent pas directement sur les paramètres adimensionés L et 

Tous les sphéroïdes obtenus se caractérisent par un moment angulaire L < 0.2 et par une vitesse angulaire < 
0.3. Les essais effectués montrent que, pour un moment angulaire L > 0.2, les corps formés sont des objets triaxiaux 
dont la position dans le plan (L, Ù^) ne coïncide pas avec les formes d'équilibres de Jacobi. On pense [29] que 
les ellipsoïdes ainsi obtenus pourraient, sur des échelles de temps plus longues, connaître un ajustement de forme 
ainsi qu'une diminution de leur vitesse de rotation, évoluant ainsi progressivement vers des formes d'équilibres 
de Jacobi. Au-delà de L ~ 0.2, l'effondrement gravitationnel entre en compétition avec la rotation rapide du 
nuage. Ainsi, selon la forme du nuage et le profil de vitesse utilisé, l'effondrement aura lieu uniquement pour la 



Par formation directe, on entend la formation - via rcaccumulation - de 2 astéroïdes rcaccumulcs en orbite l'un autour de l'autre. 
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partie centrale du nuage. Un nuage très allongé peut alors se fragmenter en sous-ensembles pour donner naissance 
à plusieurs corps pouvant, au gré des interactions, fusionner ou former des systèmes binaires ou multiples. Ce 
scénario semble se vérifier dans les simulations : les satellites obtenus sont tous issus de l'effondrement de nuages 
très aplatis, vérifiant dx/dy = d^/d^ — 16. Ce comportement, favorale à la formation de satellites, est semblable 
à celui rencontré par Michel et al. [5]. 

De plus, dans le plan {L, H^), deux corps pourvus d'un identique occuperont des positions différentes, au- 
trement dit des L différents, si et seulement si leurs formes respectives sont différentes. 

En effet, un ellipsoïde tournant autour de son axe d'inertie maximale (son petit axe), possède un moment angulaire 
L — Ilû où le moment d'inertie / s'écrit 7 = (a^ -I- b^). 

Il est important de signaler que les conditions initiales ne contiennent aucune "signature" a priori de la forme 
finale qu'on obtiendra en fin de simulation. En effet, ceci est attesté par la variété des conditions imposées qui, 
pour différents moments angulaires L, aboutissent à des objets très similaires. 

Par exemple, afin de "forcer" le corps à adopter une forme différente de la forme sphéroïdale, plusieurs simulations 
ont été réalisés en déformant sensiblement le nuage initial : sphérique à l'origine, le nuage a été aplati pour devenir 
ellipsoïdal. Les corps formés au cours de ces simulations se situent majoritairement aux environs du sommet de la 
séquence de Maclaurin. Caractérisés par un L plus élevé que précédemment (0.2 < L < 0.4) et par ces 
corps possèdent des périodes de rotation variant entre ~ 71i et ~ 91i, cohérent avec les périodes observées pour les 
astéroïdes évoluant dans la Ceinture Principale. 

D'autre part, ces simulations semblent indiquer qu'un nuage ellipsoïdal, doté d'un moment angulaire suffisant, 
est propice à la formation d'un ou plusieurs corps réaccumulés s'organisant autour d'un corps central. En effet, 
on observe dans la zone définie auparavant (0.2< L <0.4 ; >0.4) plusieurs de cas de formation de satellites 
d'astéroïdes. 

4.2. Formation de satellites 

Parmi les simulations effectuées, quatre cas ont abouti à la formation de corps secondaires (composés de plus 
d'une particule) en orbite autour d'un corps primaire. Ces corps secondaires ont tous été produits à partir de 
nuages très fortement allongés {dx — 4x10'* m; dy = dz — 2.5x10^ m) pour lesquels il a été observé une frag- 
mentation en plusieurs aggrégats au cours de l'effondrement gravitationnel. Hormis cinq cas oii le pourcentage de 
masse réaccrétée atteint 100%, toutes les simulations montrent des particules en orbite du corps réaccumulé. 

Parmi ces simulations (tab. 3) : 

HUl : formation d'un compagnon (corps secondaire) à partir d'une perte de masse du corps primaire. Le rapport 
de masse secondaire / primaire est tel que Ms/Mrp = 0.086). A titre de comparaison, on a pour l'astéroïde double 
Herminone Ms/Mrp ~ 0.09 

H9 : formation de 2 corps réaccumulés dont l'un s'échappe rapidement : absence de satellisation. Le rapport de 
masse secondaire / primaire est faible : Ms/Mrp = 0.008. 

Hll : formation de 3 corps réaccumulés. Alors que l'un des 3 est éjecté, les 2 autres se mettent en orbite l'un 
autour de l'autre. Avec un rapport Ms/Mrp = 0.859, ce couple fait fortement penser à l'astéroïde binaire Antiope 
dont les deux composantes sont comparables en taille et en masse. 

H12 : formation de 2 corps réaccumulés en orbite l'un autour de l'autre. Ici encore, avec un rapport Ms/Mrp = 
0.792, ce système binaire obtenu numériquement rappelle l'astéroïde Antiope. 



5. Discussion 

Les formes d'équilibres obtenues dans l'étude de Richardson et al. [29] occupent une bonne partie du plan 
(L, Q/^) et se situent toutes sans exception à l'intérieur d'un domaine de stabilité défini par la relation (3). Toutes 
ces formes ont également pour particularité de se situer en-deça de la limite théorique calculé par Holsapple [28] 
pour un angle de frottement 4> ~ 40° ce qui renforce encore l'idée que les astéroïdes possèdent une structure de 
type "agglomérée" , en opposition à une structure entièrement monolithique. 

De notre côté, nous observons une variété de formes possibles davantage réduite par rapport à [29]. Cela signi- 
fierait donc que, si en effet une grande variété de formes est possible (ce qui a déjà été démontré par la stabilité 
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des corps formés par [29]), seul un sous-ensemble d'entre elles peuvent être le résultat de la seule réaccumulation 
gravitationnelle. Les simulations effectuées suggèrent en effet que le processus de réaccumulation ne suffit pas à 
reproduire la grande variété de formes réellement observée dans la population d'astéroïdes. Cependant, nos simu- 
lations ont permis de montrer que des satellites sont susceptibles d'être créés lors de ces phases de réaccumulation 
gravitationnelle, confirmant ainsi ce qui avait déjà été observé dans Michel et al. [5] avec une utilisation différente 
du code numérique utilisé ici. Les satellites formés au cours de nos simulations résultent d'une fragmentation en 
plusieurs agrégats au cours l'effondrement du nuage. Il n'a pas été possible, en revanche, d'observer la formation 
de satellites via la fission d'un corps réaccumulé. 

En effet, lorsque le moment angulaire L devient trop élevé, la théorie des formes d'équilibre prévoit, comme 
solution stable, un objet binaire en orbite serrée, formé de deux composantes identiques. Dès lors, on pourrait 
s'attendre à observer un gros corps qui se divise, ou une accumulation aboutissant directement à deux compo- 
santes. 

On pense que le fait de n'avoir pas pu observer de fission est une conséquence de la démarche utilisée pour mener 
les simulations. Une fois le corps formé et stabilisé, à la suite de l'effondrement gravitationnel du nuage, il n'y a a 
priori aucune raison pour que son moment angulaire augmente soudainement pour entraîner la fission du corps. 
Un tel phénomène est plus susceptible d'être observé en adoptant une démarche inverse à la notre : appliquer des 
vitesses de rotation croissantes à un corps réaccumulé déjà formé. De plus, l'orbite des satellites obtenus au cours 
des simulations n'est pas définitive. Forces de marées, champ de gravité du Soleil sont en effet autant de facteurs 
extérieurs pouvant influencer l'évolution de l'orbite des satellites sur des échelles de temps beaucoup plus longues 
que la simulation en elle-même. 

Toutefois, il a été récemment rémarqué [27] que les astéroïdes les plus proches de la séquence d'équilibre fluide 
étaient en réalité des systèmes binaires. Parmi ces corps primaires, aucun n'est un ellipsoïde de Jacobi (à 3 axes), 
alors que tous se sont révélés être proches de sphéroïdes aplatis (ellipsoïdes à 2 axes ou sphéroïdes de Maclaurin) . 

Un résultat analogue se dégage de nos simulations d'effondrement gravitationnel : il semble en effet très dif- 
ficile d'obtenir des corps triaxiaux alors que, par opposition, les simulations forment presque exclusivement des 
sphéroïdes aplatis quelque soit le moment angulaire L du nuage initial. Comme évoqué en début de partie, ce 
résultat semble indiquer que les astéroïdes triaxiaux que l'on observe ne peuvent être formés directement par 
réaccumulation gravitationnelle. En effet, tout astéroïde subit au cours du temps de nombreux impacts mineurs, 
certes insuffisants pour être détruit, mais dont on pense qu'ils seraient suffisamment énergétiques pour modifier 
sensiblement leur forme ainsi que leur vitesse de rotation, leur procurant ainsi une position différente dans le plan 

Les sphéroïdes aplatis constitueraient ainsi une catégorie de formes dites "primordiales" qui, au gré d'éventuels 
impacts, pourraient se remodeler et évoluer vers des formes triaxiales, tels certains astéroïdes que l'on observe. 
Il convient toutefois de noter que ce mécanisme de modification des formes par impacts successifs n'a pas encore 
été étudié. Il sera donc intéressant d'analyser les formes observées par photométrie lors des prochaines grandes 
campagnes d'observations (Gaia, Pan-STARRS), car selon des études récentes [33] une partie des objets d'un 
diamètre supérieur à 20 km pourrait être des objets dits "primordiaux", c'est-à-dire n'ayant pas subi d'impacts 
très énergétiques depuis leur formation. 

Plus généralement, s'il s'avère que les formes biaxiales dominent dans la formation par réaccumulation gravita- 
tionnelle, on peut alors s'attendre à ce que le pourcentage de sphéroïdes soit plus important parmi les objets dont 
la fréquence d'impact est plus petite. Une telle vérification sera possible lorsque les données seront en nombre suf- 
fisant. D'autre part, il est tout aussi intéressant de noter que les corps binaires formés par simulations se trouvent 
dans la même région du plan {L,Ù^) que les corps binaires observés. Une telle observation suggère en effet que les 
corps binaires observés pourraient s'être formés par un processus similaire à la réaccumulation gravitationnelle. 



6. Conclusion 

Au cours de ce travail, nous avons montré que, même si une grande variété de formes peut exister pour des 
agrégats auto-gravitants, il est possible que seulement une partie d'entre elles soient accessibles par réaccumulation 
gravitationnelle. La formation directe de corps avec un moment d'inertie élevé semble ne pas être favorisée, voire 
même être totalement inacessible. 

Les limitations du modèle numérique (particules sphériques et de mêmes tailles) ne peuvent toutefois être négligées 
et rendent nécessaires des tests ultérieurs, effectués par exemple avec d'autres choix sur la distribution de tailles. 



8 



de vitesses initiales, du nombre de particules, etc. De plus, il faut remarquer que tous les éventuels phénomènes 
agissant sur des échelles de temps plus longues (comme les forces de marées) ne sont pas pris en compte dans les 
simulations en raison des contraintes liées à leurs durées. 

Nous avons également confirmé que la présence de petits satellites et d'astéroïdes binaires (composantes de tailles 
similaires) était une issue possible du processus de réaccumulation, essentiellement pour des valeurs élevées de 
moment angulaire. 
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Tab. 1 

Paramètres communs à toutes les simulations. 6 : pas; nSteps : nombre de pas; e± : restitution perpendiculaire; ey : restitution 
parallèle; Mtotal '■ masse du nuage initial; pp : densité des particules; Ntot '■ nombre de particules constituant le nuage. 
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Tab. 2 

Simulations conduisant à la formation d'un corps unique (sans satellites). L est le moment angulaire (unité : 
Mq UA^ (années/27r)~^) ; Nj-p/Ntot la fraction de particules composant le corps réaccumulé; b/a et c/a les rapports d'axes; 
e l'aplatissement ; p^p la densité et P la période de rotation du corps formé. Les simulation Del, De2, De3 ont été réalisée avec pp 
— 2g cm~^, Pp — 3 g cm~^, pp — 4 g cm~^ 
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Tab. 3 

Simulations conduisant à la formation d'un unique corps accompagné d'un ou plusieurs satellites. L est le moment angulaire (unité : 
Mq .UA^ .(années/27r) ""*"), N^p/Ntot la fraction de particules composant le corps réaccumulc, b/a et c/a les rapports d'axes, e 
l'aplatissement, p^p la densité et P la période de rotation du corps forme. 

Simulation Taille Nuage L N^p/Ntot b/a c/a e p^p P Ns Ah-^ 

{profil de vitesses) {x,y,z X 10"^ m) (x 10~^^ {kg.m^'"*) 
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FiG. 1. Courbes de Maelaurin et Jacobi (traits pleins). La partie de la courbe de Maclaurin située à droite du point de bifur- 
cation (Maclaurin - Jacobi) correspond à des états d'équilibres instables. Les états d'équilibres d'un fluide compressible (n est 
l'indice polytropiquc) sont également réprésentés (tirets). L'aplatissement théorique des planètes (ro7ids), supposées homogènes et 
incompressibles, est supérieur aux valeurs observées (croix). D'après [31]. 
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FiG. 2. De gauche à droite, de haut en bas : effondrement gravitationnel d'un nuage sphcrique (simulation DIO, cf. tableau 2) 
conduisant à la formation d'un corps ré-accumulé sans satellite. Chaque image correspond respectivement à i — 0.01 jours, t — 0.17 
jours, t — 0.23 jours, t — 0.29 jours, t — 0.47 jours et i — 23.25 jours (fin de simulation). L'objet final présente un petit nombre de 
particules-satellites. Sa forme est celle d'un sphéroïde aplati dont l'axe de rotation est pratiquement perpendiculaire au plan de la 
figure. 
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FiG. 3. Rapports d'axes obtenus par observation photomctriquc (à gauche) et des corps formes au cours des simulations numériques 
(à droite). Triangles violets : résultats de l'effondrement gravitationnel d'un nuage sphérique (d^ — dy — d^), carrés bleus 
turquoises : même chose pour un nuage aplati (dx — dy), cercles verts : même chose pour un nuage aplati {dy — dz)-, étoiles 
rouges : effondrement d'un nuage sphérique mais pour pp ^ Ig cm~^, disques bleus : cas de formation de satellites. La séqenee de 
Maclaurin et de Jacobi sont représentées, respectivement, en vert et en rouge. 
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FiG. 4. Les différents projections de l'espace des paramètres de forme et rotation, défini par a j c^h j ■X'- Triangles violets : 
résultats de l'effondrement gravitationnel d'un nuage sphcrique {d^ — dy — ds), carrés bleus turquoise : même chose pour un nuage 
aplati (dx — dy), cercles verts : même chose pour un nuage aplati {dy — rf^:), étoiles rouges : effondrement d'un nuage sphérique 
mais avec pinit > Ig cm^^, disques bleus ; cas de formation de satellites. La séqence de Maclaurin et de Jacobi sont représentées, 
respectivement, en bleu et en rouge. 
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